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1. Introduction 

It has been known for several decades that Dirac's bra-ket formalism is mathematically 
justified not by the Hilbert space alone, but by the rigged Hilbert space (RHS). This 
is the reason why there is an increasing number of Quantum Mechanics textbooks that 
already include the rigged Hilbert space as part of their contents (see, for example, 
Refs. jTj-jH])- Despite the importance of the RHS, there is still a lack of simple examples 
for which the corresponding RHS is constructed in a didactical manner. Even worse, 
there is no pedagogical discussion on the physical significance of the RHS. In this paper, 
we use the one- dimensional (ID) rectangular barrier potential to introduce the RHS at 
the graduate student level. As well, we discuss the physical significance of each of the 
ingredients that form the RHS. The construction of the RHS of such a simple model 
will unambiguously show that the RHS is needed at the most basic level of Quantum 
Mechanics. 

The present paper is complemented by a previous publication, Ref. [TUj . to which we 
shall refer the reader interested in a detailed mathematical account on the construction 
of the RHS of the ID rectangular barrier. For a general background on the Hilbert 
and the rigged Hilbert space methods, the reader may consult Ref. [TT] and references 
therein. 

Dirac's bra-ket formalism was introduced by Dirac in his classic monograph |12j . 
Since its inception, Dirac's abstract algebraic model of bras and kets (from the bracket 
notation for the inner product) proved to be of great calculational value, although there 
were serious difficulties in finding a mathematical justification for the actual calculations 
within the Hilbert space, as Dirac and von Neumann themselves state in their 
books [H|. As part of his bra-ket formalism, Dirac introduced the so-called Dirac delta 
function, a formal entity without a counterpart in the classical theory of functions. It was 
L. Schwartz who gave a precise meaning to the Dirac delta function as a functional over 
a space of test functions ^3] . This led to the development of a new branch of functional 
analysis, the theory of distributions. By combining von Neumann's Hilbert space with 
the theory of distributions, I. Gelfand and collaborators introduced the RHS [TBJIH]. R 
was already clear to the creators of the RHS that their formulation was the mathematical 
support of Dirac's bra-ket formalism [IBj- The RHS made its first appearance in the 
Physics literature in the 1960s [13120112115 when some physicists also realized that the 
RHS provides a rigorous mathematical rephrasing of all of the aspects of Dirac's bra- 
ket formalism. Nowadays, there is a growing consensus that the RHS, rather than the 
Hilbert space alone, is the natural mathematical setting of Quantum Mechanics [22]. 

A note on semantics. The word "rigged" in rigged Hilbert space has a nautical 
connotation, such as the phrase "fully rigged ship;" it has nothing to do with any 
unsavory practice such as "fixing" or predetermining a result . The phrase "rigged Hilbert 
space" is a direct translation of the phrase "osnashchyonnoe Hilbertovo prostranstvo" 
from the original Russian. A more faithful translation would be "equipped Hilbert 
space." Indeed, the rigged Hilbert space is just the Hilbert space equipped with 
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distribution theory — in Quantum Mechanics, to rig a Hilbert space means simply to 
equip that Hilbert space with distribution theory. Thus, the RHS is not a replacement 
but an enlargement of the Hilbert space. 

The RHS is neither an extension nor an interpretation of the physical principles 
of Quantum Mechanics, but rather the most natural, concise and logic language to 
formulate Quantum Mechanics. The RHS is simply a mathematical tool to extract 
and process the information contained in observables that have continuous spectrum. 
Observables with discrete spectrum and a finite number of eigenvectors (e.g., spin) do 
not need the RHS. For such observables, the Hilbert space is sufficient. Actually, as 
we shall explain, in general only unbounded observables with continuous spectrum need 
the RHS. 

The usefulness of the RHS is not simply restricted to accounting for Dirac's bra- 
ket formalism. The RHS has also proved to be a very useful research tool in the 
quantum theory of scattering and decay (see Ref. jTTj and references therein), and 
in the construction of generalized spectral decompositions of chaotic maps [23112] • In 
fact, it seems that the RHS is the natural language to deal with problems that involve 
continuous and resonance spectra. 

Loosely speaking, a rigged Hilbert space (also called a Gelfand triplet) is a triad of 
spaces 

3>cHc$ x (1.1) 

such that 7i is a Hilbert space, $ is a dense subspace of Ti and <fr x is the space 
of antilinear functionals over $ Mathematically, $ is the space of test functions, 
and <fr x is the space of distributions. The space <fr x is called the antidual space of 
Associated with the RHS (jl.lj) . there is always another RHS, 

(1.2) 

where is called the dual space of $ and contains the linear functionals over |26j . 

The basic reason why we need the spaces <J?' and <fr x is that the bras and kets 
associated with the elements in the continuous spectrum of an observable belong, 
respectively, to <£' and <fr x rather than to 7i. The basic reason reason why we need 
the space <fr is that unbounded operators are not defined on the whole of Ti but only 
on dense subdomains of 7i that are not invariant under the action of the observables. 
Such non-invariance makes expectation values, uncertainties and commutation relations 
not well defined on the whole of TC. The space $ is the largest subspace of the Hilbert 
space on which such expectation values, uncertainties and commutation relations are 
well defined. 

The original formulation of the RHS ^01 HZ] does not provide a systematic 
procedure to construct the RHS generated by the Hamiltonian of the Schrodinger 
equation, since the space <fr is assumed to be given beforehand. Such systematic 
procedure is important because, after all, claiming that the RHS is the natural setting 
for Quantum Mechanics is about the same as claiming that, when the Hamiltonian 
has continuous spectrum, the natural setting for the solutions of the Schrodinger 
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equation is the RHS rather than just the Hilbert space. The task of developing a 
systematic procedure to construct the RHS generated by the Schrodinger equation was 
undertaken in Ref. The method proposed in Ref. [H], which was partly based 

on Refs. [T[JJ 1201 EE] , has been applied to two simple three-dimensional potentials, see 
Refs. |27| 128]. to the three-dimensional free Hamiltonian, see Ref. [29J, and to the ID 
rectangular barrier potential, see Ref. [TU]. In this paper, we present the method of 
Ref. P3] in a didactical manner. 

The organization of the paper is as follows. In Sec. we outline the major reasons 
why the RHS provides the mathematical setting for Quantum Mechanics. In Sec. we 
recall the basics of the ID rectangular potential model. Section 0] provides the RHS of 
this model. In Sec. we discuss the physical meaning of each of the ingredients that 
form the RHS. In Sec. we discuss the relation of the Hilbert space spectral measures 
with the bras and kets, as well as the limitations of our method to construct RHSs. 
Finally, Sec. [7] contains the conclusions to the paper. 

2. Motivating the rigged Hilbert space 

The linear superposition principle and the probabilistic interpretation of Quantum 
Mechanics are two major guiding principles in our understanding of the microscopic 
world. These two principles suggest that the space of states be a linear space (which 
accounts for the superposition principle) endowed with a scalar product (which is used 
to calculate probability amplitudes). A linear space endowed with a scalar product is 
called a Hilbert space and is usually denoted by TC (SOj- 

In Quantum Mechanics, observable quantities are represented by linear, self- 
adjoint operators acting on 7i. The eigenvalues of an operator represent the possible 
values of the measurement of the corresponding observable. These eigenvalues, which 
mathematically correspond to the spectrum of the operator, can be discrete (as the 
energies of a particle in a box), continuous (as the energies of a free, unconstrained 
particle), or a combination of discrete and continuous (as the energies of the Hydrogen 
atom) . 

When the spectrum of an observable A is discrete and A is bounded [HI], then A 
is defined on the whole of 7i and the eigenvectors of A belong to 7i. In this case, A 
can be essentially seen as a matrix. This means that, as far as discrete spectrum is 
concerned, there is no need to extend 7i. However, quantum mechanical observables are 
in general unbounded [31] and their spectrum has in general a continuous part. In order 
to deal with continuous spectrum, textbooks usually follow Dirac's bra-ket formalism, 
which is a heuristic generalization of the linear algebra of Hermitian matrices used for 
discrete spectrum. As we shall see, the mathematical methods of the Hilbert space are 
not sufficient to make sense of the prescriptions of Dirac's formalism, the reason for 
which we shall extend the Hilbert space to the rigged Hilbert space. 

For pedagogical reasons, we recall the essentials of the linear algebra of Hermitian 
matrices before proceeding with Dirac's formalism. 
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2.1. Hermitian matrices 

If the measurement of an observable A (e.g., spin) yields a discrete, finite number N 
of results a n , n = 1, 2, . . . , N, then A is realized by a Hermitian matrix on a Hilbert 
space 7i of dimension N. Since 7^ is an iV-dimensional linear space, there are N linearly 
independent vectors {e n }^ =1 that form an orthonormal basis system for 7i. We denote 
these basis vectors e n also by \e n ). The scalar products of the elements of the basis 
system are written in one of the following ways: 

' Cm = (fin? Cm) = \Cfi,|Cm} &nm j "^i HI 1) 2, . . . , A'" , (2.1) 

where 5 nm is the Kronecker delta. As the basis system for the space 7i, it is always 
possible to choose the eigenvectors of A. Therefore, one can choose basis vectors e n ETC 
which also fulfill 

Ae n = a n e n . (2.2) 

Since A is Hermitian, the eigenvalues a n are real. The eigenvectors e n are often labeled 
by their eigenvalues a n and denoted by 

e„ = \a n ) , (2.3) 

and they are represented by column vectors. For each column eigenvector e n = \a n ), 
there also exists a row eigenvector e n = (a n \ that is a left eigenvector of A, 

e n A = a n e n . (2.4) 

Thus, when A is a Hermitian matrix acting on an A^-dimensional Hilbert space 7i, for 
each eigenvalue a n of A there exist a right (i.e., column) eigenvector of A 

A\a n ) = a n \a n ) , n = 1, 2, . . . , N , (2.5) 
and also a left (i.e., row) eigenvector of A 

(a n \A = a n (a n \ , ra = 1, 2, . . . , iV, (2.6) 
such that these row and column eigenvectors are orthonormal, 

(a n \a m )=S nm , n,m = 1,2, . . . , N, (2.7) 
and such that every vector (p GTi. can be written as 

N 

¥ = 22 \ a n)(a n \f) ■ (2.8) 

n=l 

Equation ()2.8|) is called the eigenvector expansion of ip with respect to the eigenvectors 
of A. The complex numbers (a n \<p) are the components of the vector p> with respect to 
the basis of eigenvectors of A. Physically, (a n \(p) represents the probability amplitude 
of obtaining the value a n in the measurement of the observable A on the state ip. By 
acting on both sides of Eq. (|2.8J) with A, and recalling Eq. (|2.5|) . we obtain that 

N 

Aip = y^ j a n \a n ){a n \ip) . (2.9) 

71=1 
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2.2. Dirac's bra-ket formalism 

Dirac's formalism is an elegant, heuristic generalization of the algebra of finite 
dimensional matrices to the continuous-spectrum, infinite-dimensional case. Four of 
the most important features of Dirac's formalism are: 

(i) To each element of the spectrum of an observable A, there correspond a left and a 
right eigenvector (for the moment, we assume that the spectrum is non-degenerate). 
If discrete eigenvalues are denoted by a n and continuous eigenvalues by a, then the 
corresponding right eigenvectors, which are denoted by the kets \a n ) and \a), satisfy 

A\a n ) = a n \a n ) , (2.10a) 

A|o) = a|a), (2.106) 

and the corresponding left eigenvectors, which are denoted by the bras (a n \ and (a\, 
satisfy 

(a n \A = a n (a n \ , (2.11a) 

(a\A = a(a\. (2.116) 

The bras (a\ generalize the notion of row eigenvectors, whereas the kets \a) 
generalize the notion of column eigenvectors. 

(ii) In analogy to Eq. ()2.8j) . the eigenbras and eigenkets of an observable form a complete 
basis, that is, any wave function (p can be expanded in the so-called Dirac basis 
expansion: 

<f = ^2\a n )(a n \(f) + / da\a)(a\(p) . (2.12) 

In addition, the bras and kets furnish a resolution of the identity, 

|a„)(a n | + J da|a)(a|, (2.13) 

and, in a generalization of Eq. (|2.9|) . the action of A can be written as 

A = a n \a n )(a n \ + / daa\a)(a\. (2-14) 

n J 

(iii) The bras and kets are normalized according to the following rule: 

(a n \a m ) = 5 nm , (2.15a) 

(a\a!) — 8(a — a') , (2.156) 

where 5 nm is the Kronecker delta and 5 (a — a') is the Dirac delta. The Dirac 
delta normalization generalizes the orthonormality (|2.7j) of the eigenvectors of a 
Hermitian matrix. 

(iv) Like in the case of two finite-dimensional matrices, all algebraic operations such as 
the commutator of two observables A and B, 

[A, B] = AB — BA , (2.16) 

are always well defined. 



E 
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2. 3. The need of the rigged Hilbert space 

In Quantum Mechanics, observables are usually given by differential operators. In the 
Hilbert space framework, the formal prescription of an observable leads to the definition 
of a linear operator as follows: One has to find first the Hilbert space 7i, then one sees 
on what elements of TC the action of the observable makes sense, and finally one checks 
whether the action of the observable remains in TC. For example, the position observable 
Q of a ID particle is given by 

Qf(x)=xf(x). (2.17) 

The Hilbert space of a ID particle is given by the collection of square integrable 
functions, 

/oo 
dx|/(x)| 2 <oo}, (2.18) 
-oo 

and the action of Q, although in principle well defined on every element of L 2 , remains 
in L 2 only for the elements of the following subspace: 

/oo 
dx\xf(x)\ 2 < oo}. (2.19) 
■oo 

The space T>(Q) is the domain of the position operator. Domain ()2.19|) is not the whole 
of L 2 , since the function g(x) = l/(x + i) belongs to L 2 but not to T>(Q); as well, 
Q is an unbounded operator, because \\Qg\\ = oo; as well, QV{Q) is not included in 
T>(Q), since h(x) = \j{x 2 + 1) belongs to V{Q) but Qh does not belong to V{Q). The 
denseness and the non-invariance of the domains of unbounded operators create much 
trouble in the Hilbert space framework, because one has always to be careful whether 
formal operations are valid. For example, Q 2 = QQ is not defined on the whole of L 2 , 
not even on the whole of T>(Q), but only on those square integrable functions such that 
x 2 f G L 2 . Also, the expectation value of the measurement of Q in the state (p, 

(<p,Q<p), (2.20) 

is not finite for every ip G L 2 , but only when tp G T>(Q). Similarly, the uncertainty of 
the measurement of Q in (p, 



A^Q = ^{p,Q 2 p)-{p,Qp) 2 , (2.21) 

is not defined on the whole of L 2 . 

On the other hand, if we denote the momentum observable by 

Pf(x) = -ih±f(x) , (2.22) 
dx 

then the product of P and Q, PQ, is not defined everywhere in the Hilbert space, but 
only on those square integrable functions for which the quantity 

PQf(x) = -^xf(x) = -ih (f(x) + xf'(x)) (2.23) 

makes sense and is square integrable. Obviously, PQf makes sense only when / is 
differentiate, and PQf remains in L 2 only when /, /' and xf are also in L 2 ; thus, 
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PQ is not denned everywhere in L 2 but only on those square integrable functions that 
satisfy the aforementioned conditions. Similar domain concerns arise in calculating the 
commutator of P with Q. 

As in the case of the position operator, the domain T>(A) of an unbounded operator 
A does not coincide with the whole of TC but is just a dense subspace of TC [25J; 
also, in general T>(A) does not remain invariant under the action of A, that is, AT>(A) 
is not included in T>(A). Such non-invariance makes expectation values, 

(<p,A<p), (2.24) 

uncertainties, 

A V A = y/(<p,A*ip)-(<p,A<p)*, (2.25) 

and algebraic operations such as commutation relations not well defined on the whole 
of the Hilbert space TC [34 . Thus, when the position, momentum and energy operators 
Q, P, H are unbounded, it is natural to seek a subspace <fr of TC on which all of these 
physical quantities can be calculated and yield meaningful, finite values. Because the 
reason why these quantities may not be well defined is that the domains of Q, P and 
H are not invariant under the action of these operators, the subspace $ must be such 
that it remains invariant under the actions of Q, P and H. This is why we take as <3? 
the intersection of the domains of all the powers of Q, P and H [THj: 

oo 

$= p) V{A n B m ). (2.26) 

n,m=0 
A,B=Q,P,H 

This space is known as the maximal invariant subspace of the algebra generated by Q, P 
and H, because it is the largest subdomain of the Hilbert space that remains invariant 
under the action of any power of Q, P or H, 

A*c$, A = Q,P,H. (2.27) 

On all physical quantities such as expectation values and uncertainties can be 
associated well-defined, finite values, and algebraic operations such as the commutation 
relation ()2.16|) are well defined. In addition, the elements of are represented by 
smooth, continuous functions that have a definitive value at each point, in contrast 
to the elements of TC, which are represented by classes of functions which can vary 
arbitrarily on sets of zero Lebesgue measure. 

Not only there are compelling reasons to shrink the Hilbert space TC to but, 
as we are going to explain now, there are also reasons to enlarge TC to the spaces 
<fr x and <&' of Eqs. and (jl.2j) . When the spectrum of A has a continuous part, 

prescriptions (J2.11frj) and (|2.10foj) associate a bra (a\ and a ket \a) to each element a of the 
continuous spectrum of A. Obviously, the bras (a\ and kets \a) are not in the Hilbert 
space [HE], and therefore we need two linear spaces larger than the Hilbert space to 
accommodate them. It turns out that the bras and kets acquire mathematical meaning 
as distributions. More specifically, the bras (a\ are linear functionals over the space 
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and the kets \a) are antilinear functionals over the space That is, (a\ G <&' and 
| a) e # x . 

In this way, the Gelfand triplets of Eqs. (jl.lj) and (jl.2j) arise in a natural way. 
The Hilbert space TC arises from the requirement that the wave functions be square 
normalizable. Aside from providing mathematical concepts such as self-adjointness or 
unitarity, the Hilbert space plays a very important physical role, namely 7i selects 
the scalar product that is used to calculate probability amplitudes. The subspace <fr 
contains those square integrable functions that should be considered as physical, because 
any expectation value, any uncertainty and any algebraic operation can be calculated 
for its elements, whereas this is not possible for the rest of the elements of the Hilbert 
space. The dual space and the antidual space 3> x contain respectively the bras and 
the kets associated with the continuous spectrum of the observables. These bras and 
kets can be used to expand any ip 6 $ as in Eq. (j2.12|) . Thus, the rigged Hilbert space, 
rather than the Hilbert space alone, can accommodate prescriptions (j2.10aj) - (j2.16j) of 
Dirac's formalism. 

It should be clear that the rigged Hilbert space is just a combination of the Hilbert 
space with distribution theory. This combination enables us to deal with singular objects 
such as bras, kets, or Dirac's delta function, something that is impossible if we only use 
the Hilbert space. 

Even though it is apparent that the rigged Hilbert space should be an essential 
part of the mathematical methods for Quantum Mechanics, one may still wonder if the 
rigged Hilbert space is a helpful tool in teaching Quantum Mechanics, or rather is a 
technical nuance. Because basic quantum mechanical operators such as P and Q are 
in general unbounded operators with continuous spectrum [SEj, and because this kind 
of operators necessitates the rigged Hilbert space, it seems pertinent to introduce the 
rigged Hilbert space in graduate courses on Quantum Mechanics. 

From a pedagogical standpoint, however, this section's introduction to the rigged 
Hilbert space is not sufficient. In the classroom, new concepts are better introduced 
by way of a simple, exactly solvable example. This is why we shall construct the RHS 
of the ID rectangular barrier system. We note that this system does not have bound 
states, and therefore in what follows we shall not deal with discrete spectrum. 

2.4- Representations 

In working out specific examples, the prescriptions of Dirac's formalism have to be 
written in a particular representation. Thus, before constructing the RHS of the ID 
rectangular barrier, it is convenient to recall some of the basics of representations. 

In Quantum Mechanics, the most common of all representations is the position 
representation, sometimes called the x-representation. In the x-representation, the 
position operator Q acts as multiplication by x. Since the spectrum of Q is (— oo, oo), 
the x-representation of the Hilbert space 7i is given by the space L 2 . In this paper, we 
shall mainly work in the position representation. 
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In general, given an observable B, the ^-representation is that in which the operator 
B acts as multiplication by 6, where the 6's denote the eigenvalues of B. If we denote 
the spectrum of B by Sp(B), then the 6-representation of the Hilbert space 7i is given 
by the space L 2 (Sp(£>), db), which is the space of square integrable functions f{b) with 
b running over Sp(-B). In the 6-representation, the restrictions to purely continuous 
spectrum of prescriptions ()2.10a|) - (j2.13p become 

(b\A\a) = a(b\a) , (2.28a) 
(a\A\b) = a(a\b) , (2.286) 

(b\<p) = / da(b\a)(a\ip) , (2.28c) 



5(b - b') = (b\b') = J da (b\a) (a\b') . (2.28d) 

The "scalar product" (b\a) is obtained from Eq. (|2.28aj) as the solution of a differential 
eigenequation in the 6-representation. The (b\a) can also be seen as transition elements 
from the a- to the 6-representation. Mathematically, the (b\a) are to be treated as 
distributions, and therefore they often appear as kernels of integrals. In this paper, we 
shall encounter a few of these "scalar products" such as (x\p), (x\x') and (x\E )i, r - 

3. Example: The one-dimensional rectangular barrier potential 

The example we consider in this paper is supposed to represent a spinless particle moving 
in one dimension and impinging on a rectangular barrier. The observables relevant to 
this system are the position Q, the momentum P, and the Hamiltonian H. In the 
position representation, Q and P are respectively realized by the differential operators 
(I2.17J) and (|2.22|) . whereas H is realized by 

/ h 2 d 2 \ 
Hf(x)=(- — — + V(x))f(x), (3.1) 



2m dx 2 



where 




— oo < x < a 

V(x) = { V Q a<x<b (3.2) 
b < x < oo 

is the ID rectangular barrier potential. Formally, these observables satisfy the following 
commutation relations: 

[Q,P]=iM, (3.3a) 
[H t Q] = ~P, (3.36) 

[H,P] = ih^. (3.3c) 
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Since our particle can move in the full real line, the Hilbert space on which the 
differential operators (ETTTj) . (l2~22[) and (JSUJ) should act is L 2 of Eq. (I2~TH[) . The 

corresponding scalar product is 



(f,g)= / dxf(x)g(x), f,geL 2 } (3.4) 



V{Q) = 


{feL 2 \ 


xfeL 2 } , 


V{P) = 


{/ e ^ 2 1 


/ G AC, P/GL 2 } , 


V(H) = 


{f £ L 2 


/ G AC 2 , tf/GZ 2 } 



where /(x) denotes the complex conjugate of f{x). 

The differential operators (|2.17|) . (|2.22|) and (|3.1|1 induce three linear operators on 
the Hilbert space L 2 . These operators are unbounded [TOj, and therefore they cannot 
be defined on the whole of L 2 , but only on the following sub domains of L 2 [ID) : 

(3.5a) 
(3.55) 
(3.5c) 

where, essentially, AC is the space of functions whose derivative exists, and AC 2 is the 
space of functions whose second derivative exists (see Ref. (TOJ for more details). On 
these domains, the operators Q, P and H are self-adjoint [10 . 

In our example, the eigenvalues (i.e., the spectrum) and the eigenfunctions of the 
observables are provided by the Sturm-Liouville theory. Mathematically, the eigenvalues 
and eigenfunctions of operators extend the notions of eigenvalues and eigenvectors of a 
matrix to the infinite-dimensional case. The Sturm-Liouville theory tells us that these 
operators have the following spectra |1U) : 

Sp(Q) = (— oo, oo) , (3.6a) 

Sp(P) = (-oo, oo) , (3.65) 

Sp(#) = [0, oo) . (3.6c) 

These spectra coincide with those we would expect on physical grounds. We expect the 
possible measurements of Q to be the full real line, because the particle can in principle 
reach any point of the real line. We also expect the possible measurements of P to be 
the full real line, since the momentum of the particle is not restricted in magnitude or 
direction. The possible measurements of H have the same range as that of the kinetic 
energy, because the potential does not have any wells of negative energy, and therefore 
we expect the spectrum of H to be the positive real line. 

To obtain the eigenfunction corresponding to each eigenvalue, we have to solve the 
eigenvalue equation ([2.106)1 for each observable. Since we are working in the position 
representation, we have to write Eq. ([2.106)1 in the position representation for each 
observable: 

(x\Q\x') = x'(x\x') , (3.7a) 
(x\P\p) = p(x\p) , (3.76) 
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(x\H\E) = E(x\E) . (3.7c) 

By recalling Eqs. (J2D3J), (PUZj) and (jSHJ), we can write Eqs. (jSZZ^-JSZZS) as 

a^xja/) = , (3.8a) 

— ih-^-(x\p) = p(x\p) , (3.86) 

/r <h + ^(x)^ (x|£) = E(x\E) . (3.8c) 



2m dx 2 



For each position a;', Eq. ()3.8aj) yields the corresponding eigenfunction of Q as a delta 
function, 

(x\x') = 6(x - x') . (3.9) 

For each momentum p, Eq. (|3.86|) yields the corresponding eigenfunction of P as a plane 

wave, 

pipx/h 

( ^7Sr (3 ' 10) 

For each energy E, Eq. ()3.8c|) yields the following two linearly independent 
eigenfunctions [10J: 

1 ,„ f T(k)e~ [kx -oo < x < a 

/ ffj \ 1/2 J v ' 

( x |£+) r= ^__j x<^ A r (A;)e iKa; + B r {k)e~ iKX a<x<b (3.11a) 
' R r (k)e ihx + e~ ihx b < x < oo , 



gifcz _|_ RJ^q ikx _ QO < ^ < a 

(x|£+)i= — — I x<^ A x (k)e 1KX + BAk)e- iKX a<x<b (3.116) 
\2-nkn 1 J \ m/ , N 

T(k)e lkx b < x < oo , 



where 



and where the coefficients that appear in Eqs. (|3.11a|l - (|3.116jl can be easily found by 
the standard matching conditions at the discontinuities of the potential JU|- Thus, in 
contrast to the spectra of Q and P, the spectrum of H is doubly degenerate. 

Physically, the eigenfunction (x\E + ) r represents a particle of energy E that impinges 
on the barrier from the right (hence the subscript r) and gets reflected to the right with 
probability amplitude R r {k) and transmitted to the left with probability amplitude T(k), 
see Fig. EJi. The eigenfunction (x\E + )\ represents a particle of energy E that impinges 
on the barrier from the left (hence the subscript 1) and gets reflected to the left with 
probability amplitude R\(k) and transmitted to the right with probability amplitude 
T(k), see Fig. Eh . 
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Note that, instead of (|3.11aj) - (|3.11fij) . we could choose another pair of linearly 
independent solutions of Eq. ()3.8cj) as follows [TU] : 

, ,„ f T*(k)e [kx -oo < x < a 

f 771 \ 

(x\E-) v = {■^jtfp) x \ A*(k)e- iKX + B;{k)e iKX a<x<b (3.13a) 

I R*(k)e~ ikx + e ikx b < x < oo , 



e lkx + Rf(k)e ikx — oo < x < a 

(x\Er)i = [yT^J X ) A *( k ) e ~ iKX + B*(k)e iKX a<x<b (3.135) 



m \ V 2 

T*(k)e~ ikx b<x<oo 



where the coefficients of these eigenfunctions can also be calculating by means of the 
standard matching conditions at x = a, b [TU] . The eigenfunction (x\E~) r represents 
two plane waves — one impinging on the barrier from the left with probability amplitude 
T*(k) and another impinging on the barrier from the right with probability amplitude 
R* (k) — that combine in such a way as to produce an outgoing plane wave to the right, see 
Fig. |2^. The eigenfunction (x\E~)\ represents two other planes waves — one impinging 
on the barrier from left with probability amplitude R*(k) and another impinging on the 
barrier from the right with probability amplitude T*(k) — that combine in such a way 
as to produce an outgoing wave to the left, see Fig. The eigensolutions (x\E~) Tj i 
correspond to the final condition of an outgoing plane wave propagating away from the 
barrier respectively to the right and to the left, as opposed to (x\E + ) Tt i, which correspond 
to the initial condition of a plane wave that propagates towards the barrier respectively 
from the right and from the left. 

The eigenfunctions (jUJ), (J3~T0|) . (j3~TT^ - (j3~TUj ) and (j3~T3^ - (l3~T3l are not square 
integrable, that is, they do not belong to L 2 . Mathematically speaking, this is the reason 
why they are to be dealt with as distributions (note that all of them except for the delta 
function are also proper functions). Physically speaking, they are to be interpreted in 
analogy to electromagnetic plane waves, as we shall see in Section 



4. Construction of the rigged Hilbert space 

In the previous section, we saw that the observables of our system are implemented by 
unbounded operators with continuous spectrum. We also saw that the eigenfunctions 
of the observables do not belong to L 2 . Thus, as we explained in Sec. we need to 
construct the rigged Hilbert spaces of Eqs. (jl.lj) and (jl.2j) [see Eqs. (I4.8j) and (j4.21j) 
below]. We start by constructing 

4-1. Construction of $ = «S(R-{a, b}) 

The subspace $ is given by Eq. ()2.26|) . In view of expressions ()2.17|) . ()2.22|) and (|3.1|) . 
the elements of $ must fulfill the following conditions: 

• they are infinitely differentiable, so the differentiation operation can be applied as 
many times as wished, 
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• they vanish at x = a and x = b, so differentiation is meaningful at the discontinuities 
of the potential [3*T] . 

• the action of all powers of Q, P and H remains square integrable. 
Hence, 

$ = { v e L 2 | <p e C°°(R), p {n \a) = <^ (n) (&) = , n = 0, 1, ... , 

P n Q m H l y{x) eL 2 , n, m, I = 0, 1, . . .} , (4.1) 

where C°°(R) is the collection of infinite different iable functions, and denotes the 
nth derivative of ip. From the last condition in Eq. (|4.1J) . we deduce that the elements 
of <fr satisfy the following estimates: 



(p\\ n , m ,i = JJ dx \P n Q m H l ip{x)\ 2 < oo, n, m, I = 0, 1, . . . . (4.2) 

These estimates mean that the action of any combination of any power of the observables 
remains square integrable. For this to happen, the functions <f(x) must be infinitely 
differentiate and must fall off at infinity faster than any polynomial. The estimates (|4.2jl 
induce a topology on <fr, that is, they induce a meaning of convergence of sequences, in 
the following way. A sequence {(p a } ^-converges to ip when {(p a } converges to (p with 
respect to all the estimates (|4.2j) . 

Pa-^P if \\p a ~ <p\\n,m,l ►O, 71, 171, I = 0, 1, . . . . (4.3) 

q^oo a— +oo 

Intuitively, a sequence p a converges to ip if whenever we follow the terms of the sequence, 
we get closer and closer to the limit point ip with respect to a certain sense of closeness. 
In our system, the notion of closeness is determined by the estimates || \\ ntm ,i, which 
originate from the physical requirements that led us to construct 

From Eqs. (j4.1|) and ()4.2j) . we can see that $ is very similar to the Schwartz space 
the major differences being that the derivatives of the elements of $ vanish at 
x = a,b and that $ is not only invariant under P and Q but also under H. This is why 
we shall write 

$ = S(M-{a,6}) . (4.4) 

It is always a good, though lengthy exercise to check that «S(M.-{a, b}) is indeed 
invariant under the action of the observables, 

A5(l-{fl,fc})c5(l-{a,6}), A = P,Q,H. (4.5) 

This invariance guarantees that the expectation values 

(<p, A n ip) , ip e 5(R-{a, b}) , A = P,Q,H, n = 0, 1, . . . (4.6) 

are finite, and that the commutation relations (|3.3aj) - (l3.3c|) are well defined [HE]. It 
can also be checked that P, Q and H, which are not continuous with respect the 
topology of the Hilbert space L 2 , are now continuous with respect to the topology 

of s(R-{a, b}) maim. 
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4-2. Construction o/$ x = 5 x (R-{a, b}). The Dirac kets 

The space <fr x is simply the collection of r^-continuous antilinear functionals over <E> [2t 
By combining the spaces 7i and <!>*, we obtain the RHS of our system, 

# c n c $ x , (4.7) 

which we denote in the position representation by 

<S(R-{a, b}) CL 2 cS x (R-{a, 6}) . (4.8) 

The space iS x (R-{a,6}) is meant to accommodate the eigenkets \p), \x) and lE^)^ of 
P, Q and if. In the remainder of this subsection, we construct these eigenkets explicitly 
and see that they indeed belong to S x (R-{a, b}). We shall also see that \p), \x) and 
l-E*)^ are indeed eigenvectors of the observables. 

The definition of a ket is borrowed from the theory of distributions as follows |16j . 
Given a function f(x) and a space of test functions the antilinear functional F that 
corresponds to the function f(x) is an integral operator whose kernel is precisely f(x): 



F(<p) = J dxcp(x)f(x), (4.9a) 
which in Dirac's notation becomes 

(<p\F) = [ dx(<p\x)(x\f). (4.95) 



It is important to keep in mind that, though related, the function f(x) and the functional 
F are two different things, the relation between them being that f(x) is the kernel of F 
when we write F as an integral operator. In the physics literature, the term distribution 
is usually reserved for f(x). 

Definition ()4.9aj) provides the link between the quantum mechanical formalism 
and the theory of distributions. In practical applications, what one obtains from the 
quantum mechanical formalism is the distribution f(x) (in this paper, the plane waves 
_=^ e ipzM ; the delta function S(x — x') and the eigenfunctions (xlE^i^). Once f(x) is 
given, one can use definition ()4.9aj) to generate the functional \F). Then, the theory of 
distributions can be used to obtain the properties of the functional \F), which in turn 
yield the properties of the distribution f(x). 

By using prescription (|4.9»| . we can define for each eigenvalue p the eigenket \p) 
associated with the eigenfunction (j3.10j) : 

/°° i 
dx^x)-==e^ h , (4.10a) 

which, using Dirac's notation for the integrand, becomes 

/■oo 

{(p\p)= / dx(<p\x)(x\p) . (4.106) 



Similarly, for each x, we can define the ket \x) associated with the eigenfunction ()3.9j) 
of the position operator as 



(cp\x) = / dx' <f(x')5(x - x') , (4.11a) 
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which, using Dirac's notation for the integrand, becomes 

POO 

{(f\x) = / da;' ((p\x'){x'\x) . (4.116) 



The definition of the kets l-E^)^ that correspond to the Hamiltonian's eigenfunc- 
tions (|3.11a|) - (|3.116J) and (|3.13fljl - (|3.13fo)) follows the same prescription: 



(^1^)1.= / dx^(a;)(x|E ± ) 1 , r , (4.12a) 



{ip\E^)^ = / dxiip^ixlE^. (4.126) 



that is, 



(Note that this equation defines four different kets.) One can now show that the 
definition of the kets \p), \x) and l-E^)^ makes sense, and that these kets indeed belong 
to the space of distributions S x (R-{a, b}) [TO] . 

As in the general case of Eqs. ()4.9aj) - fl4.96|) . it is important to keep in mind the 
difference between eigenfunctions and kets. For instance, (x\p) is an eigenfunction of a 
differential equation, Eq. (j3.86|) . whereas \p) is a functional, the relation between them 
being given by Eq. (j4.106|) . A similar relation holds between (x'\x) and \x), and between 
(x\E )i tT and lE^)^. It is also important to keep in mind that "scalar products" like 
(x\p), (x'\x) or (xlE^ij do not represent an actual scalar product of two functionals; 
these "scalar products" are simply solutions to differential equations. 

We now turn to the question of whether the kets \p), \x) and \E^)i tI are eigenvectors 
of the corresponding observable [see Eqs. ()4.17|) - (j4.19j) below]. Since the observables act 
in principle only on their Hilbert space domains, and since the kets lie outside the Hilbert 
space, we need to extend the definition of the observables from <fr into <fr x , in order to 
specify how the observables act on the kets. The theory of distributions provides us 
with a precise prescription of how an observable acts on $ x , and therefore of how it 
acts on the kets, as follows [TSJ. The action of a self-adjoint operator A on a functional 
\F) e <fr x is defined as 

(<p\A\F) = (A<p\F) , for all <p in $ . (4.13) 

Note that this definition extends the Hilbert space definition of a self-adjoint operator, 

(f,Ag) = (Af,g), (4.14) 

which is valid only when / and g belong to the domain of A. In turn, Eq. ()4.13|) can 
be used to define the notion of eigenket of an observable: A functional \a) in <fr x is an 
eigenket of A with eigenvalue a if 

(ip\A\a) = (A<f\a) = a(ip\a) , for all (p in $ . (4.15) 

When the "left sandwiching" of this equation with the elements of $ is understood and 
therefore omitted, we shall simply write 

A\a) = a\a), (4.16) 
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which is just Dirac's eigenket equation (|2.106|) . Thus, Dirac's eigenket equation acquires 
a precise meaning through Eq. ()4.15|) . in the sense that it has to be understood as "left 
sandwiched" with the wave functions ip of 

By using definition ()4.15|h one can show that \p), \x) and l-E^)^ are indeed 
eigenvectors of P, Q and H, respectively fTU| : 

P\p)=p\p), peR, (4.17) 
Q\x)=x\x), xeR, (4.18) 
HIE^ = ElE^, EG[0,oo). (4.19) 

4-3. Construction of = S'(R-{a,b}) . The Dirac bras 

In complete analogy with the construction of the Dirac kets, we construct in this 
subsection the Dirac bras (p\, (x\ and i, r ( ± -E| of P, Q and H. Mathematically, the 
Dirac bras are distributions that belong to the space which is the space of linear 
functionals over $ The corresponding RHS is 

$ C H C , (4.20) 

which we denote in the position representation by 

S(R-{a, b}) C L 2 C <S'(R-{a, b}) . (4.21) 

Likewise the definition of a ket, the definition of a bra is borrowed from the theory 
of distributions jTHj. Given a function f(x) and a space of test functions the linear 
functional F generated by the function f(x) is an integral operator whose kernel is the 
complex conjugate of f(x): 



F(<p) = j dx<p(x)f(x), (4.22a) 
which in Dirac's notation becomes 

(F\p)= [ dx(f\x)(x\ V ). (4.226) 



Note that this definition is very similar to that of a linear functional, Eq. (|4.9ojl . except 
that the complex conjugation affects f(x) rather than <f(x), which makes F linear rather 
than antilinear. Likewise the antilinear case (|4.9ojl . it is important to keep in mind that, 
though related, the function f(x) and the functional F are two different objects, the 
relation between them being that f(x) is the kernel of F when we write F as an integral 
operator. 

By using prescription (|4.22aJ) . we can now define for each eigenvalue p the eigenbra 
(p\ associated with the eigenfunction ()3.10|) : 

r°° i 

(p\if}= / dxip(x)^=e- [px/h } (4.23a) 



V2vrn 

which, using Dirac's notation for the integrand, becomes 

POO 

(p\<p)= / dx(p\x)(x\tp) . (4.236) 
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Comparison with Eq. (|4.1Uaj) shows that the action of (p\ is the complex conjugate of 
the action of \p), 



(P\V) = (V\P) , (4-24) 

and that 

(p\x)=Wi=-^=e-^/ h . (4.25) 
The bra (x\ is defined as 

/oo 

dx' (f(x')5(x — x 1 ) , (4.26a) 
■oo 

which, using Dirac's notation for the integrand, becomes 

(x\f) = / dx' (x\x')(x'\(p) . (4.266) 



Comparison with Eq. ()4.11aj) shows that the action of (x\ is complex conjugated to the 
action of \x), 



(x\<p) = (<p\x) , (4.27) 

and that 

(x\x') = (x'\x) = 5(x — x) . (4.28) 
Analogously, the eigenbras of the Hamiltonian are defined as 

/oo 
dx ip(x) u^Elx) , (4.29a) 
■oo 

that is, 



/oo 
dx^EWixltp), (4.296) 
-oo 



where 



i ir f^|ar> = (x\E% T . (4.30) 

(Note that in Eq. ()4.29aj) we have defined four different bras.) Comparison of Eq. (j4.29oj) 
with Eq. (j4.12q|) shows that the actions of the bras i, r ( ± -E| are the complex conjugates 
of the actions of the kets \E )\/- 

l , v ( ± E\cp)=J^) ltI . (4.31) 

Now, by using the RHS mathematics, one can show that the definitions of (p|, (x\ and 
i ir ( ± i?| make sense and that (p\, (x\ and i, r ( ± -E| belong to 5'(M-{a,6}) [THj . 

Our next task is to see that the bras we just defined are left eigenvectors of the 
corresponding observable [see Eqs. ()4.35j) - (j4.37j) below]. For this purpose, we need to 
specify how the observables act on the bras, that is, how they act on the dual space 
^'(^-{a, 6}). We shall do so in analogy to the definition of their action on the kets, by 
means of the theory of distributions jTH] . The action to the left of a self-adjoint operator 
A on a linear functional (F\ G 3>' is defined as 

(F\A\<p) = (F\A<p) , for all in $ . (4.32) 
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Likewise definition (|4.13|) . this definition generalizes Eq. (J4.14j) . In turn, Eq. (|4.32j) can 
be used to define the notion of eigenbra of an observable: A functional (a\ in is an 
eigenbra of A with eigenvalue a if 

(a\A\(p) = (a\A(p) = a(a\(p) , for all y? in $ . (4.33) 

When the "right sandwiching" of this equation with the elements of # is understood 
and therefore omitted, we shall simply write 

(a\A = a(a\ , (4.34) 

which is just Dirac's eigenbra equation (j2.11fc|) . Thus, Dirac's eigenbra equation acquires 
a precise meaning through Eq. ()4.33|) . in the sense that it has to be understood as "right 
sandwiched" with the wave functions <p of 

By using definition ()4.33j) . one can show that (p\, (x\ and i, r ( ± -E'| are indeed left 
eigenvectors of P, Q and H, respectively [TUj : 

(p\P = p(p\, pGR, (4.35) 

(x\Q = x(x\, xeR, (4.36) 

i^ElH = E , Ee[0,oo). (4.37) 

It is worthwhile noting that, in accordance with Dirac's formalism, there is a one- 
to-one correspondence between bras and kets [22] ; that is, given an observable A, to 
each element a in the spectrum of A there correspond a bra (a\ that is a left eigenvector 
of A and also a ket \a) that is a right eigenvector of A. The bra (a\ belongs to 
whereas the ket \a) belongs to 3> x . 

4-4- The Dirac basis expansions 

A crucial ingredient of Dirac's formalism is that the bras and kets of an observable form 

a complete basis system, see Eqs. (|2.12j) and (j2.13|) . When applied to P, Q and H, 
Eq. ijinTty yields 

dp\p)(p\ = I, (4.38) 

oo 

dx'\x')(x'\ = J, (4.39) 



poo poo 

/ dE\E ± ) u ( ± E\+ dE\E ± ) TT ( ± E\= I , (4.40) 
Jo Jo 

In the present subsection, we derive various Dirac basis expansions for the algebra of 
the ID rectangular barrier potential. We will do so by formally sandwiching Eqs. ()4.38jl - 
1)4.40)1 in between different vectors. 

If we sandwich Eqs. (J4.38)) - (J4.40)) in between (x\ and (p, we obtain 

POO 

(x\<p) = / dp(x\p){p\ip) , (4.41) 
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POO 

(x\ip) = / dx' (x\x')(x'\<p) , (4.42) 



-oo 
oo 



(x\<p)= / dE(x\E ± ) u ( ± E\ V ) + / dE(x\E ± ) vv { ± E\ V ). (4.43) 
Jo Jo 

Equations (j4.41j) - (j4.43j) can be rigorously proved by way of the RHS [10J. In proving 
these equations, we give meaning to Eqs. (j4.38j) - (j4.40J) . which are just formal equations: 
Equations fj4.HSjl - fj4.40jl have always to be understood as part of a "sandwich." Note 
that Eqs. f|4.41j) - ()4.43|) are not valid for every element of the Hilbert space but only 
for those <p that belong to 5(IR-{a, &}), because the action of the bras and kets is well 
defined only on 5(IR-{a,6}) jlHj. Thus, the RHS, rather than just the Hilbert space, 
fully justifies the Dirac basis expansions. Physically, the Dirac basis expansions provide 
the means to visualize wave packet formation out of a continuous linear superposition 
of bras and kets. 

We can obtain similar expansions to Eqs. (j4.41j) - (j4.43j) by sandwiching Eqs. (j4.38jl - 
f!4.40j) in between other vectors. For example, sandwiching Eq. f|4.39|) in between (p\ 
and <p yields [TU] 



(p\ip) = dx (p\x)(x\(p) , (4.44) 

J — oo 

and sandwiching Eq. (j4.39j) in between i, r ( ± -E , | and ip yields [ J!?] 

POO 

^Elip) = / dx^EWixltp). (4.45) 



It is worthwhile noting the parallel between the Dirac basis expansions and the Fourier 
expansions (|4.41j) and (|4.44jl [10J. This parallel will be used in Sec. El to physically 
interpret the Dirac bras and kets. 

We can also sandwich Eqs. fj4.38jl - fj4.40jl in between two elements ip and ip of 
S(R-{a,b}), and obtain PU 



M)= / dp(ip\p)(p\ij} } (4.46) 

J — oo 

POO 

M)= / dx(ip\x)(x\^), (4.47) 



roo poo 

{ip,ip)= dE{ V \E ± ) u { ± E\^)+ dE^E^^l^E^). (4.48) 
Jo Jo 

Equations fj4.46jl - fj4.48p allow us to calculate the overlap of two wave functions ip and ip 
by way of the action of the bras and kets on those wave functions. 

The last aspect of Dirac's formalism we need to implement is prescription f|2.14jl . 
which expresses the action of an observable A in terms of the action of its bras and kets. 
When applied to P, Q and H, prescription (j2.14jl yields 

dpp\p)(p\, (4.49) 
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POO 

Q = dxx\x)(x\ , (4.50) 



oo 
oo 



H= / dE E\E )u( E\ + / dEE\E ± ) vv ( ± E\ . (4.51) 
Jo Jo 

Needless to say, these equations are formal expressions that acquire meaning when 
properly sandwiched. For example, sandwiching them in between (x\ and (p yields 



(x\P(p) = / dpp(x\p) (p\(f) , (4-52) 

J — oo 

POO 

(x\Q(f) = / dx' x' (x\x'){x'\<p) , (4.53) 



oo 
oo 



(x\H<p) = / d££(x|£ ± ) 11 ( ± £|^>+ / dE E(x\E ) IT { E\(p) , (4.54) 
Jo 

and sandwiching them in between two elements (p and t/> of 5(lR-{a, b}) yields [TU] 

/•oo 

(<p,Pij})= dpp{y\p){p\ip) , (4.55) 



-oo 
oo 



((p,Qip) = / drra;(<^|x)(x|'0) , (4.56) 



/•oo /*oo 

(ip,H*l;)= / d££(^|£ ± > 11 ( ± £|V;>+ / dE ^(^| J B ± ) r • (4.57) 
Jo Jo 

Note that, in particular, the operational definition of an observable — according to which 
an observable is simply an operator whose eigenvectors form a complete basis such that 
Eqs. (I2.12J) . ()2.13|) and ()2.14j) hold, see for example Ref. [JT] — acquires meaning within 
the RHS. 

The sandwiches we have made so far always involved at least a wave function <p of 
5(lR-{a, b}). When the sandwiches do not involve elements of 5(M-{a, b}) at all, we 
obtain expressions that are simply formal. These formal expressions are often useful 
though, because they help us understand the meaning of concepts such as the delta 
normalization or the "matrix elements" of an operator. Let us start with the meaning 
of the delta normalization. When we sandwich Eq. ()4.39|) in between (p'\ and \p), we 
get 

/■oo 

dx (p'\x)(x\p) = (p'\p) . (4.58) 

This equation is a formal expression that is to be understood in a distributional sense, 
that is, both sides must appear smeared out by a smooth function ip{p) = (p\<p) in an 
integral over p: 

poo POO POO 

dpip(p) / dx (p' \x) (x\p) = / dpip(p)(p'\p) . (4.59) 



The RHS in Quantum Mechanics 22 
The left-hand side of Eq. (|4.59|) can be written as 

/»oo /»oo poo 

dx(p'\x) / dp(f(p)(x\p) = / dx{p'\x) / dp (x\p) (p\(f) 



oo 
oo 



dx (p'\x) (x\<p) 
= <f{p') (4.60) 

Plugging Eq. ()4~fiU)l into Eq. (|Q5|t leads to 

By recalling the definition of the delta function, we see that Eq. ()4.6H) leads to 

(p'\p) = 5(p-p>), (4.62) 

and to 

PCX) 

dx (p'\x)(x\p) = 5{p — p) . (4.63) 

By using Eq. 1)4.25)1 . we can write Eq. 1)4.63)1 in a well-known form: 
1 f°° 

— / dx e^ p - p ' )x/h = 8(p - p') . (4.64) 
2tt/i ,/_ 00 

This formal equation is interpreted by saying that the bras and kets of the momentum 
operator are delta normalized. That the energy bras and kets are also delta normalized 
can be seen in a similar, though slightly more involved way [2*H] : 

a^E'lE^fs = 8(E - E') 5 Qf3 , (4.65a) 

/oo 
dx a ( ± E'|x)(x|E ± )^ = 5{E-E')5 af3 , (4.656) 
-oo 

where a, (3 stand for the labels 1, r that respectively denote left and right incidence. 
The derivation of expressions involving the Dirac delta function such as Eqs. (J4.62)) . 
(14.64)1 or (|4.65ajl - (|4. 656)1 shows that these formal expressions must be understood in a 
distributional sense, that is, as kernels of integrals that include the wave functions <p of 
S(R-{a,b}), like in Eq. (jOHjl . 

In a similar way, we can also understand the meaning of the "matrix elements" of 
the observables in a particular representation, e.g.: 

(x\Q\x'} =x'5(x-x'), (4.66) 

(x\P\x) = -ih4- Six - x') , (4.67) 
dx 

(x|#|x'> = (—-fx + V(x)) 5(x - x') . (4.68) 
\ 2m dx z J 

Equations ()4. 66)1 - ()4. 68)1 can be obtained by formally inserting Eq. (j4.39j) into respectively 
Eq. (J2T7J), ([I22D and Q . 
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It is illuminating to realize that the expressions (|4.66|) - (j4.68|) generalize the matrix 



representation of an observable A in a finite-dimensional Hilbert space. If a\, 
are the eigenvalues of A, then, in the basis { |ai) , . . . , |ajv)}, A is represented as 

/ ai •■■ \ 
a 2 ■■■ 



A 



\ 







a N J 



(4.69) 



which in Dirac's notation reads as 

(aj|^4|oj) = aiSij . (4.70) 

Clearly, expressions (|4.fifi|) - (|4.fi8j) are the infinite-dimensional extension of expres- 
sion (pETOj) . 



5. Physical meaning of the Dirac bras and kets 

The bras and kets associated with eigenvalues in the continuous spectrum are not 
normalizable. Hence, the standard probabilistic interpretation does not apply to 
them straightforwardly. In this section, we are going to generalize the probabilistic 
interpretation of normalizable states to the non-normalizable bras and kets. As well, in 
order to gain further insight into the physical meaning of bras and kets, we shall present 
the analogy between classical plane waves and the bras and kets. 

In Quantum Mechanics, the scalar product of the Hilbert space is employed to 
calculate probability amplitudes. In our example, the Hilbert space is L 2 , and the 
corresponding scalar product is given by Eq. ()3.4j) . That an eigenvalue of an observable 
A lies in the discrete or in the continuous part of the spectrum is determined by this 
scalar product. An eigenvalue a n belongs to the discrete part of the spectrum when its 
corresponding eigenfunction f n (x) = (x\a n ) is square normalizable: 



/oo 
AX \fn(x)\ 2 < OO 
■oo 



(5.1) 



An eigenvalue a belongs to the continuous part of the spectrum when its corresponding 
eigenfunction f a (x) = (x\a) is not square normalizable: 



if ay fa] 



dx\f a (x)f 



oo . 



(5.2) 



In the latter case, one has to use the theory of distributions to "normalize" these states, 
e.g., delta function normalization: 



if ay fa' 



dxf a (x)f a >(x) =5(a-a'). 



(5.3) 



This Dirac delta normalization generalizes the Kronecker delta normalization of 
"discrete" states: 

POO 

(fnyfn')= / dx f n (x)f n >(x) = 5 nn > . (5.4) 
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Because they are square integrable, the "discrete" eigenvectors f n (x) = (x\a n ) can be 
interpreted in the usual way as probability amplitudes. But because they are not square 
integrable, the "continuous" eigenvectors f a (x) = (x\a) must be interpreted as "kernels" 
of probability amplitudes, in the sense that when we multiply (x\a) by (ip\x) and then 
integrate, we obtain the density of probability amplitude (<p\a): 

/oo 
dx ((p\x) (x\a) . (5.5) 
-oo 

Thus, in particular, (x\p), (x\x') and (x\E ± )i jI represent "kernels" of probability 
amplitudes. 

Another way to interpret the bras and kets is in analogy to the plane waves 
of classical optics and classical electromagnetism. Plane waves e lkx represent 
monochromatic light pulses of wave number k and frequency (in vacuum) w = kc. 
Monochromatic light pulses are impossible to prepare experimentally; all that can be 
prepared are light pulses <p(k) that have some wave-number spread. The corresponding 
pulse in the position representation, (p(x), can be "Fourier decomposed" in terms of the 
monochromatic plane waves as 

cp(x) = ^= [ dke ikx V (k) , (5.6a) 



/2tt 

which in Dirac's notation becomes 

(x\<p) = J dk(x\k)(k\<p) . (5.66) 

Thus, physically preparable pulses can be expanded in a Fourier integral by the 
unpreparable plane waves, the weights of the expansion being <p(k). When ip(k) is highly 
peaked around a particular wave number k$, then the pulse can in general be represented 
for all practical purposes by a monochromatic plane wave e lk ° x . Also, in finding out how 
a light pulse behaves under given conditions (e.g., reflection and refraction at a plane 
interface between two different media), we only have to find out how plane waves behave 
and, after that, by means of the Fourier expansion (j5.6flj) . we know how the light pulse 
<p(x) behaves. Because obtaining the behavior of plane waves is somewhat easy, it is 
advantageous to use them to obtain the behavior of the whole pulse |42j . 

The quantum mechanical bras and kets can be interpreted in analogy to the classical 
plane waves. The eigenfunction (x\p) = e ipx '^/\/zifh represents a particle of sharp 
momentum p; the eigenfunction (x\x') = 8{x — x 1 ) represents a particle sharply localized 
at x'\ the monoenergetic eigenfunction (x\E )i T represents a particle with well-defined 
energy E (and with additional boundary conditions determined by the labels ± and 1, r). 
In complete analogy to the Fourier expansion of a light pulse by classical plane waves, 
Eq. ()5.6aj) . the eigenfunctions (x\p), (x\x') and (xlE^)^ expand a wave function ip, see 
Eqs. (I4.41jl - (|4.43jl . When the wave packet (pip) is highly peaked around a particular 
momentum p , then in general the approximation (p(x) ~ e vp ° x / h / v2nh holds for all 
practical purposes; when the wave packet (fix) is highly peaked around a particular 
position x , then in general the approximation (p(x) ~ Six — x ) holds for all practical 
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purposes; and when (p(E) is highly peaked around a particular energy E , then in general 
the approximation (p(x) ~ (x\E^)i )Z holds for all practical purposes (up to the boundary 
conditions determined by the labels ± and 1, r). Thus, although in principle (x\p), (x\x') 
and (a;|i? ± )i ir are impossible to prepare, in many practical situations they can give good 
approximations when the wave packet is well peaked around some particular values po, 
Xo, E of the momentum, position and energy. Also, in finding out how a wave function 
behaves under given conditions (e.g., reflection and transmission off a potential barrier), 
all we have to find out is how the bras and kets behave and, after that, by means of 
the Dirac basis expansions, we know how the wave function ip(x) behaves. Because 
obtaining the behavior of the bras and kets is somewhat easy, it is advantageous to use 
them to obtain the behavior of the whole wave function [33] . 

From the above discussion, it should be clear that there is a close analogy between 
classical Fourier methods and Dirac's formalism. In fact, one can say that Dirac's 
formalism is the extension of Fourier methods to Quantum Mechanics: Classical 
monochromatic plane waves correspond to the Dirac bras and kets; the light pulses 
correspond to the wave functions (p; the classical Fourier expansion corresponds to the 
Dirac basis expansions; the classical Fourier expansion provides the means to form light 
pulses out of a continuous linear superposition of monochromatic plane waves, and the 
Dirac basis expansions provide the means to form wave functions out of a continuous 
linear superposition of bras and kets; the classical uncertainty principle of Fourier 
Optics corresponds to the quantum uncertainty generated by the non-commutativity 
of two observables [44J. However, although this analogy is very close from a formal 
point of view, there is a crucial difference from a conceptual point of view. To wit, 
whereas in the classical domain the solutions of the wave equations represent a physical 
wave, in Quantum Mechanics the solutions of the equations do not represent a physical 
object, but rather a probability amplitude — In Quantum Mechanics what is "waving" 
is probability. 



6. Further considerations 



In Quantum Mechanics, the main objective is to obtain the probability of measuring 
an observable A in a state (p. Within the Hilbert space setting, such probability can be 
obtained by means of the spectral measures E a of A (see, for example, Ref. |Hj). These 
spectral measures satisfy 

dE a (6.1) 

S P (A) 

and 

A = [ a dE a (6.2) 

JSp(A) 

Comparison of these equations with Eqs. ()2.13|) and (j2.14|) yields 

dE a = \a)(a\ da . (6.3) 
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Thus, the RHS is able to "factor out" the Hilbert space spectral measures in terms of 
the bras and kets [15]. For the position, momentum and energy observables, Eq. ()6.3|) 
reads as 

dE x = \x)(x\ dx , (6.4) 

dE p = \p)(p\ dp, (6.5) 

dE E = {E^u^EldE +\E ± ) rr ( ± E\dE . (6.6) 

Although the spectral measures dE a associated with a given self-adjoint operator A are 
unique, the factorization in terms of bras and kets is not. For example, as we can see 
from Eq. (|6.6jh the spectral measures of our Hamiltonian can be written in terms of the 
basis {\E + )i tT } or the basis {\E~)\ tT }. From a physical point of view, those two basis are 
very different. As we saw in Sec. the basis {|i? + )i >r } represents the initial condition 
of an incoming particle, whereas the basis {l-E^)!^} represents the final condition of an 
outgoing particle. However, the spectral measures of the Hilbert space are insensitive to 
such difference, in contrast to the RHS, which can differentiate both cases. Therefore, 
when computing probability amplitudes, the RHS gives more precise information on 
how those probabilities are physically produced than the Hilbert space. 

In this paper, we have restricted our discussion to the simple, straightforward 
algebra of the ID rectangular barrier. But, what about more complicated potentials? 
In general, the situation is not as easy. First, the theory of rigged Hilbert spaces as 
constructed by Gelfand and collaborators is based on the assumption that the space <E> 
has a property called nuclearity [TEJ[T7]. However, it is not clear that one can always 
find a nuclear space $ that remains invariant under the action of the observables. 
Nevertheless, Roberts has shown that such $ exists when the potential is infinitely 
often different iable except for a closed set of zero Lebesgue measure [TH]. Second, the 
problem of constructing the RHS becomes more involved when the observable A is 
not cyclic [TSj. And third, solving the eigenvalue equation of an arbitrary self-adjoint 
operator is rarely as easy as in our example. 

7. Summary and conclusions 

We have used the ID rectangular barrier model to see that, when the spectra of the 
observables have a continuous part, the natural setting for Quantum Mechanics is the 
rigged Hilbert space rather than just the Hilbert space. In particular, Dirac's bra-ket 
formalism is fully implemented by the rigged Hilbert space rather than just by the 
Hilbert space. 

We have explained the physical and mathematical meanings of each of the 
ingredients that form the rigged Hilbert space. Physically, the space 3> = «S(R-{a, b}) 
is interpreted as the space of wave functions, since its elements can be associated 
well-defined, finite physical quantities, and algebraic operations such as commutation 
relations are well defined on Mathematically, $ is the space of test functions. The 
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spaces <£' = <S'(R-{a, b}) and <& x = S x (R-{a, 6}) contain respectively the bras and 
kets associated with the eigenvalues that lie in the continuous spectrum. Physically, the 
bras and kets are interpreted as "kernels" of probability amplitudes. Mathematically, 
the bras and kets are distributions. The following table summarizes the meanings of 
each space: 



Space 


Physical Meaning 


Mathematical Meaning 




Space of wave functions tp 


Space of test functions tp 


n 


Probability amplitudes 


Hilbert space 




Space of kets a) 


Antidual space 




Space of bras (a\ 


Dual space 



We have seen that, from a physical point of view, the rigged Hilbert space does 
not entail an extension of Quantum Mechanics, whereas, from a mathematical point of 
view, the rigged Hilbert space is an extension of the Hilbert space. Mathematically, the 
rigged Hilbert space arises when we equip the Hilbert space with distribution theory. 
Such equipment enables us to cope with singular objects such as bras and kets. 

We have also seen that formal expressions involving bras and kets must be 
understood as "sandwiched" by wave functions tp. Such "sandwiching" by tp's is 
what controls the singular behavior of bras and kets. This is why mathematically the 
sandwiching by </?'s is so important and must always be implicitly assumed. In practice, 
we can freely apply the formal manipulations of Dirac's formalism with confidence, since 
such formal manipulations are justified by the rigged Hilbert space. 

We hope that this paper can serve as a pedagogical, enticing introduction to the 
rigged Hilbert space. 
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Figure 1. Schematic representation of the eigenfunctions (x\E + } r , Fig. la, and 
(x\E+) u Fig. lb. 
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Figure 2. Schematic representation of the eigenfunctions (x\E ) r , Fig. 2a, and 
(x\E~) h Fig. 2b. 



